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Abstract 

Using a unified approach of optical-mechanical analogy in a semiclassical 
formula, we evaluate the effect of Chern-Simons modified gravity on the quan- 
tum phase shift of de Broglie waves in neutron interferometry. The phase shift 
calculated here reveals, in a single equation, a combination of effects coming 
from Newtonian gravity, inertial forces, Schwarzschild and Chern-Simons mod- 
ified gravity. However the last two effects, though new, turn out to be too tiny 
to be observed, and hence only of academic interest at present. The approxima- 
tions, wherever used, as well as the drawbacks of the non-dynamical approach 
are clearly indicated. 

PACS numbers(s): 95.85.Ry, 12.15.Ff, 04.50. +h 



I. Introduction 



Investigations at a semiclassical level of gravitational effects on quantum 
interference are expected to be useful for the development of a theory of quan- 
tum gravity [1]. Although the full theory is yet to appear, outcomes of early 
seminal experiments by Colella, Overhauser and Werner (COW) [2] and an 
improved later experiment by Werner, Staudenmann and Colella (WSC) [3], 
measuring respectively the effects of Newtonian gravity and Coriolis force on 
the phase shift of neutrons, continue to provide essential resources in this di- 
rection. The effect of Coriolis force was first suggested by Page [4]. Several 
other related works have also contributed useful information. For instance, a 



1 



scmiclassical treatment reveals that quantum uncertainty in the source variables 
induces uncertainties in the metric components of gravity in a specific manner 
[5]. Another important phenomenon is the neutrino flavor oscillation induced 
by gravity [6-9]. A different kind of theoretical approach in the calculations 
of the phase of quantum particles, neutrinos included, has yielded a very in- 
teresting result: The Dirac spin 1/2 particle has the exact covariant Stodolsky 
[10] phase S/h = (1/h) J p^dx^ in a static gravity field [11]. Some of these 
investigations could provide appropriate theoretical backgrounds for astrophys- 
ical applications, e.g., in the atmospheric neutrino experiments [12], or in the 
observations involving 7-ray bursts [13]. 

On the other hand, the developments in matterwave interferomctry have 
shown a greater promise (over photon interferometry) in the measurements of 
the phase shift due to Earth's gravity and axial rotation. In the absence of 
external forces, the interference with itself of the de Broglie wave associated 
with an ensemble of particles allows us to predict, via Huygcn's principle, the 
motion of the wave. The presence of an external field modifies the motion 
causing a shift in the interference fringes. A classic example is the Aharonov- 
Bohm (AB) fringe shift which provides information as to how the motion of 
electrons are modified in the presence of a magnetic potential. However, the 
situation becomes more complicated when one introduces the effects of gravity 
and rotation on matterwave interferometry. 

An intuitive and elegant special relativistic treatment of the quantum phase 
shift for thermal neutrons has been proposed by Anandan [14]. However, the ef- 
fects of gravity and rotation are still considered as separate components. On the 
other hand, a more comprehensive analysis involving gravity and inertial forces 
would require the use of complete geometrical framework of curved spacetime. 
The special relativistic effects should follow only as a limiting case. Anandan's 
semiclassical treatment paves the way for the application of simpler "Newto- 
nian" form of equations of optical-mechanical analogy in curved spacetime [15- 
19]. These equations provide a substitute for geometrical framework though 
only for calculational purposes. Using them, the corrections to the phase shift 
had been worked out [20] in certain spacetimes: A "Coriolis" force leading to 
a gravitational analog of the quantum AB effect appeared naturally from the 
application of a "Lorentz force" form of the geodesic equation in Kerr and Kerr- 
Sen [21] spacetime. (Optical analogs are artifacts designated by and they 
can be easily translated into actual mechanical quantities, see below.) 

A different genre of theory that leads to rotating solution different from that 
of Kerr is Chern-Simons (hereafter CS) modified gravity. A reasonable summary 
of the timeline and study of CS modified gravity is in order. The CS modified 
gravity descends from all versions of string theory upon 4-dimcnsional com- 
pactification due to the Green-Schwarz anomaly-cancelling mechanism [21,22]. 
The first group to study CS modified gravity as a string inspired correction to 
general relativity consisted of Alexander, Pcshkin and Sheikh- Jabbari [23] as 
well as Alexander and Gates [24]. The initial non-dynamical formulation of CS 
modified gravity was presented by Jackiw and Pi [25] who also showed that the 
Schwarzschild metric as well as a special type of gravitational wave is a solution 
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of CS modified gravity. Guarrera and Hariton [26] showed, using Jackiw-Pi CS 
coupling field, that Reissner-Nordstrom is also a solution of CS gravity. The 
first spinning black hole solutions were found by Alexander and Yunes [27, 28]. 
These solutions represent spinning black holes in the far field limit but are not 
axisymmetric. Smith, Erickcek, Caldwell, and Kamionkowski [29] evaluated 
constraints on the CS parameter space from current satellite experiments such 
as LAGEOS and Gravity Probe-B. Their approach can be regarded as an exten- 
sion of the Jackiw-Pi formulation because the action contains a scalar field that 
is time varying but spatially homogeneous. Grumiller and Yunes [30] performed 
a detailed study of exact solutions that could represent spinning black holes with 
arbitrary CS coupling fields. This is the most complete search for black hole- like 
solutions currently available in CS modified gravity for a non-dynamical field. 
Alexander and Yunes [31] studied the coupling of fcrmions to CS modified grav- 
ity. Yunes and Sopuerta [32] studied perturbations of Schwarzschild black holes 
in CS modified gravity. 

Konno, Matsuyama, Asano and Tanda [33] (hereafter KMAT) studied a so- 
lution for a slowly rotating black hole in CS modified gravity derived previously 
by Konno, Matsuyama and Tanda [34] (hereafter KMT). The KMT slow rotat- 
ing solution corresponds to a (special) spatial coupling field. The solution differs 
from all others in that it is axisymmetric, among other things. However, the 
solution and the non-dynamical approach have several unphysical features. The 
solution is based on a purely ad hoc non-dynamical scalar, which plagues the 
well-posedness of the initial value problem [32]. Moreover, the solution is non- 
unique and does not carry over to the dynamical sector because the scalar field 
has infinite energy, which in turn prevents a self-consistent perturbation of the 
general relativity solutions [35]. Despite these drawbacks, the Alexander- Yunes 
far field solution and the KMT slow rotating solution are useful for understand- 
ing the CS effects on physical observables. KMAT [33] showed that their solution 
remarkably explains galactic flat rotation curves without the need to invoke hy- 
pothetical dark matter. They also calculated the effect of frame dragging on 
spin precession, which could provide a way to observationally test CS gravity. 

In this paper, we shall consider the KMT solution [33,34] to study the effects 
of CS modified gravity on the phase shift in neutron interferometry using, in 
addition, the equations of optical-mechanical analogy for rotating sources [36]. 
The final formula for phase shift shows a combination of different effects. We 
also recover, as a limiting case, the observed zeroth order effects of Newtonian 
gravity and Coriolis force. As shown by KMAT, one of the CS constants appear 
as constant tangential velocity in the galactic halo. Thus one might treat its 
appearance in the phase shift as a non-local effect similar to the effect of, say, 
cosmological constant in local experiments. Our calculation shows that the 
effect of CS constants, which represents dragging of inertial frames, appear only 
in the second order of slow rotation parameter and thus too tiny. 

The paper is organized as follows: To be reasonably informative, we briefly 
state in Sec. II the CS action together with the KMT solution and discuss sta- 
bility of circular orbits, which is of crucial importance. The involved approxi- 
mations on the geodesic equation are clearly mentioned. Sec. Ill describes the 
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essence of Anandan's argument which provides a core clement of the ensuing 
analysis. To familiarize the readers at large with what we refer to as optical me- 
chanical analogy in curved spacetime, we briefly review its contents in Sec. IV. 
We calculate the exact phase shift in CS gravity in Sec.V and extract the correc- 
tions in Sec. VI. In Sec. VII, we summarize the results of the paper together with 
a brief review of the difficulties of the KMT solution and the non-dynamical 
approach. We shall take G = cq = 1, unless specifically displayed. 

II. The CS action, KMT solution and stability 

(a) CS modified gravity follows from the action [33,34] 

I=^j^{y=gR-\v.K^, (1) 

where = <9 M i? is an external 4- vector, the so-called the embedding coordinate. 
The CS topological current is given by 



(2) 



where e ^"/ 3 the Levi-Civita tensor density of weight one. Variation of the action 
with respect to the metric leads to the CS field equation 

where = R^ v — \g^ v R is the Einstein tensor, is the usual energy- 
momentum tensor, and C^ v is the Cotton tensor given by 



= \ 
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v a (e a » af3 V a R" p + e^^VaR^ + (V a v T ) ^R T ^ V +*R ruat ") 

(4) 

In the above, the dual Riemann tensor is defined by *R T J iv = ^e^ va P R T aa p and 
v Tcr = \7 a v T = dudrd — Y^dx'd is a symmetric tensor. The Bianchi identity 
VjjG' 11 ' = and the equation of motion V = together imposes the 
condition [33,34] 

= V^C"" = ^=v»*R%^R^ x . (5) 

(b) KMAT [33] consider perturbation around the Schwarzschild solution of 
mass M in the small expansion parameter e(= J /Mr) where J is the angular 
momentum. Then, for a slowly rotating star, the KMT solution [34] in the 
case of a spacelike vector = d^-d = d^z = <9 M (rcos#) = (0, cos 6, — r sin#, 0) 
becomes 

ds 2 = - (l - ^ dt 2 + (l " ^) 1 dr2 + r2 ( d ° 2 + si ° 2 W*) -lr 2 w{r)dtd<j>, 

(6) 
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the ^-dependence of the (£0)-component is restricted by Eq. (5), and 



u)(r) = -4 ( 1 ] + -4[r 2 - 2Mr - 4M 2 + 4M(r - 2M) ln(r - 2M)1 (7) 



where C\ and C2 are constants (hereafter called CS constants) characteriz- 
ing slow rotation i.e., Ci,C2 ~ O(e). From the metric, it is clear that ui <~ 
(length) -1 , which implies that C\ ~ (length) and C2 is truly dimensionless. 

We should clearly mention that we are considering only the non-dynamical 
formulation of CS modified gravity, where the CS scalar field in the KMT so- 
lution satisfies the evolution equation only to leading order in the spin [35] 
(See Ref. [29] for the dynamical formulation). However, it is well to remember 
that the non-dynamical theory have several physical drawbacks as mentioned 
in the introduction. The solution (6) is based on the choice of a particular 
non-dynamical scalar defined by $ — rcos#/A , where Ao is a constant. On 
the other hand, in the non-dynamical framework, there are no well-motivated 
physical reasons for particular choices of the scalar field except for simplifying 
the equations. Furthermore, for a given choice of CS scalar the Pontryagin con- 
straint significantly restricts the class of allowed solutions, even to the point 
where the non-dynamical theory may be over-constrained and lack a well-posed 
initial value problem [32] . The hope to imbed the KMT solution in the dynam- 
ical sector is dashed because the stress energy associated with •& = rcos#/Ao 
is infinite. As shown by Yunes and Pretorius [35], even a more general choice 
of the KMT scalar does not heal the disease. Therefore the KMT metric is 
not a self-consistent solution to the dynamical field equations. Yunes and Pre- 
torius [35] have found a new stationary axisymmetric solution for a different 
non-dynamical scalar. Although the canonical $ is not allowed by this family 
of solutions, it does bypass the problem with the stress tensor and is compati- 
ble with the dynamical framework. The existence of new independent solution 
shows that a solution in the non-dynamical theory such as that of KMT can 
not be a unique black hole solution. 

In the CS modified theory, as well as in general relativity, the standard 
geodesic description holds if and only if the back-reaction is neglected in the 
equations of motion. Such an approximation is valid, for example, in the 
extreme-mass ratio limit. It is not valid, for instance, in the strong field of 
equal mass binaries. Only under the assumption that the effect of test particle's 
self-force on its own trajectory be neglected, the geodesic calculations lead to 
the transverse circular velocity v 1 ^ in the equatorial plane (9 = n/2) as 



where u' = dw/dr. The first term clearly comes from the Schwarzschild metric 
while the second term becomes a constant at large r, i.e., 



v' 




+ ruj(r) + yw'(r) 



(8) 




(9) 
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This result remarkably explains the observed flat rotation curves in the galactic 
halo. Essentially C2 is responsible for frame dragging that decays as r~ 4 . 

(c) It is however not a priori evident if the circular orbits in the spacetime 
are stable. On the other hand, stability of circular orbits in the galactic halo is 
a crucial requirement if any model of the halo has to be physically viable. We 
have carried out the necessary calculations with the potential 

nr) = -(l-^) (l + - 2£?Mr), (10) 

where the constants E and L are respectively the conserved relativistic energy 
and angular momentum per unit rest mass of the test particle in a circular 
orbit. Putting in the relevant expressions [33] for E and L in \ r —R, we 
find the following results: The expression for V" = \r=R at any arbitrary 
but fixed radius r = R depends only on C\ and C2, and that \ t =r< 
only for values of these constants very small relative to the length scale M, for 
instance C\,C<i G [0.001,0.0002] (see Fig.l for details), which of course is in 
perfect accordance with the slow rotation approximation. So we conclude that 
the orbits are indeed stable for very small values of Ci, C2. Otherwise, the 
orbits are unstable. 

III. Special relativistic phase shift: Anandan's formula 

Before we develop the equations for CS gravity, it is useful to have a brief 
preview of the special relativistic phase shift because it will ultimately provide 
the limiting case. (All special relativistic quantities will be denoted by tilde). 
Using intuitive arguments, Anandan [14] proposed a general formula for the 
quantum phase shift based on a correspondence between the shift and the clas- 
sical special relativistic equations of motion. Consider two de Broglie wavelets 
originating at A and interfering at B, one travelling along ADB and the other 
along ACB, the plane ADBC being rotated by 180° about the horizontal axis 
AO (Fig. 2) while the Earth itself is rotating about its axis with angular veloc- 
ity f2. The vertical direction is determined by the resultant of gravity and the 
centrifugal force of Earth. In the absence of any external forces, BD=BC and 
we take CD = d, and AB = I. The external fields cause a shift in the angular 
position of B by 69 ^ and thereby cause a phase shift. Let v be the velocity 
of the classical particle and 6v be the change in the velocity in the plane of 
interference. Since 69 is small, sin<5# ~ 69 ~ 4p Then the special relativistic 
phase shift A(f> SR can be calculated as [14]: 

_ Sv kddv± I pA dv± jmAdv±/dt , 

= Kd ~ = T~dTi = m?~dT = — m — ' 

where 7 = (l — w 2 /c§) , A = Id is the planar area (assumed small) enclosed 
by the two paths of interfering beams, m is the neutron mass, p{= hk) is the 
momentum in flat space, and dv±/dt is the component of acceleration in the 
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vertical direction CD. This remarkable general formula, Eq.(ll), at once gives 
rise to several individual components. 

The following cases were considered: (i) A particle with charge e moving in 
a magnetic field B. Then, jdv±/dt = eBv/m, which, when used in Eq.(ll), 
immediately yields the Aharonov-Bohm (AB) effect, viz., 

A0 AB = eBA/h. (1 2) 

(ii) A spinless particle of mass to in a gravitational field, dv±/dt = —g, where g is 
the gravitational acceleration on the surface of the Earth. Using the expression 
for the relativistic momentum j 2 = 1 +jj 2 /to 2 Cq and the de Broglie relation 
p = hk, Eq.(ll) leads to the following result: 

A0 ff = ~ g Am 2 /h 2 R - gAk/cl. (13) 

The COW experiment was accurate enough to measure the first quantum term, 
but not the second, that is, the so called relativistic term, (iii) A particle 
moving in the Coriolis force field of Earth so that, to first order in f2, dv±/dt = 

x v = 2fl n v, where tt n is the component of Earth's angular velocity Q 

normal to the apparatus. Using the Planck-Einstein law E/cq = tvjj = tticq^, 
one finds from Eq.(ll) that 

(»> 

Co 



Using the dispersion relation 



wc find from Eq.(14) that 



u 2 ~n 2 = Hp. (15) 



~ 2mAtt n hAVt n k 2 

The WSC experiment has tested the first nonrelativistic term to within a good 
accuracy, (iv) This effect comes from the coupling of particle's spin to the back- 
ground curvature. Again using Eq.(ll) together with the Mathisson-Papapetrou 
force [37], the shift comes out to be 

. , hGMAu , „, 

- -mar> (16) 

where G is the Newtonian gravitational constant, M is the gravitating mass 
and R is the distance from the center. (We shall not discuss this spin-curvature 
component in this article.) With all the above in view, we proceed to familiarize 
the readers with the salient features of our approach in the curved spacetime 
regime. 
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IV. Optical-Mechanical equations in curved spacetime 

The central idea is to introduce an optical-mechanical analogy in curved 
spacetime. The exact optical form of geodesic equations in the spherically sym- 
metric case of general relativity was originally derived in Ref. [15] . The equations 
provide an excellent tool that enables one to visualize the problems of geomet- 
rical optics as problems of classical mechanics and vice versa. The first step 
is to find out an effective optical refractive index n that is formally equivalent 
to the geometrized gravity field. This step in itself is not new. Usually the 
index equivalent to the exterior Schwarzschild field is approximately taken to 
be n w 1 + 2MG/rcQ, where M is the central gravitating mass. We derive the 
exact expression as below. 

Consider a static, spherically symmetric, but not necessarily vacuum, solu- 
tion of general relativity written in isotropic coordinates 

ds 2 = h(r)c 2 dt 2 -$- 2 {r)\dff , (17) 

where r = (x,y,z) or (r, 6, ip), and ft, $ could be the solution of Einstein's 
field equations. Many metrics of physical interest can be put into this isotropic 
form. The coordinate speed of light c(r) is determined by the condition that 
the geodesic be null {ds 2 = 0) : 

= c ^Vh, (18) 

which immediately provides an effective index of refraction for light in the grav- 
itational field given by 

„(f) _ jlj. (19) 

The concept of optical mechanical analogy can be used to recast the geodesic 
equation for both massive and massless particles into a single, exact Newtonian 
"F = ma" type of equation given by [16] 

% = V (1*4). Sift = nfrlf^f.M = * ,20) 

where m is the rest mass of the particle, Eq is the conserved total energy, V is the 
gradient operator. All the standard geodesic equations in Schwarzschild gravity 
including some new insights in cosmology follow from the above equation. This 
remarkably simple feature of the geodesic equations is brought about by the use 
of the stepping parameter A, having dimension of length, first introduced by 
Evans and Roscnquist [19]. Eqs. (20) provide an easy way to introduce quantum 
relations so that the geodesic motion of a particle can be interpreted as motion 
of de Broglie matter waves in a dispersive medium with an effective index of 
refraction N(r,X) = n(r)/Wl + (A/A c ) 2 where A c is the Compton wavelength 
given by A c = 2irh/mco and A = A<& -1 is the physical wavelength measured 



c(r) 



df 
~dt 
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in a gravity field. This interpretation allows us to extend in a straightforward 
manner the classical optical-mechanical analogy into the quantum regime. 

Alsing [36] has subsequently extended the method to broader class of metrics 
in general relativity and his equation is going to provide the basic tool in what 
follows. Consider the most general form of the metric given by 

ds 2 = goocldt 2 + 2g ai c dtdx t + gijdx l dx\ i,j — 1, 2, 3. (21) 

Define the proper time dr, proper length dl and the velocity v measured with 
respect to this proper time as 

.• dx l 



(22) 



ds 
co 



dr = —(c dt - gidx 1 ), g m = h, (23) 



v=0 



CO 



dl 2 = lt] dx l dx 1 = y- gij + ^2ij dx l dx\ v 2 = 7 y vV, (24) 

9i = -—, g i = i ij 9i = -g 0i - (25) 

#00 

The angular momentum gi and the proper velocity Vi are vectors defined in 
the 3-space characterized by the metric 7^ , which is used to raise or lower the 
indices of these 3-vectors. Now, metric (21) can be rewritten as 

ds 2 = h(c dt - g l dx i ) 2 (1 - v 2 /c 2 ) . (26) 

and the conserved energy E is given by [38] 

E = mc 2 g 0a — = a = 0,1, 2, 3. (27) 

ds ^1 - v 2 lc 2 

The variational principle for the geodesies following from Eq.(21) is given by 

5?2,*2 X2,t2 X2,t2 

S I mcods =5 j Ldt ~5 J mc^\Jh{r)\J 1 — v 2 {r,v)/c 2 j3(r 1 v)dt = 0, 

(28) 



Xl,tl 

where 



co dt pVh V ' 

From the Lagrangian L, let us find the momenta conjugate to v l . (Flat space 
quantities are designated by tilda). This is given by 

' eL + Il). (30) 



dv' V* VhJ' 
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Since the vectors gi and Vi are denned in a space with the metric 7^, we can 
identify the right hand side of the Eq.(30) as a vector in the same 3-space. Let 
us call it the momentum 3-vector 

*.Js(* + i). (31) 

It can also be verified that H — §§rv l — L = pit 1 — L = E, which is a constant 
along the trajectory of a particle as stated in Eq.(27). Hence, it is possible 

s 2 

to introduce Maupertuis principle 5 J piV l dt — 0. Assuming further that the 

£1 

spatial part of the metric could be written in an isotropic form dl = cK_e/<I>, d/e = 
5ijdx l dx° , 7^ = 5ij$~ 2 , we get 

v 2 =n 2 v 2 /f,n=^j= (32) 

and the Maupertuis variational principle yields, after lengthy algebra introduc- 
ing the parameter A, the geodesic equation in the form 

d 2 f ^±fn 2 v 2 \ dr _ x 

dA 2 = V {—) + dA X ^ X ^ ^ 

dr 

dA = n- 2 fidt, — = n 2 Vh (34) 
dJ\ 

where g = (gi). This Newtonian form of the geodesic equation, valid for both 
massless and massive particles, has been obtained under the only assumption 
that the spatial part of the metric could be written in an isotropic form. On 
eliminating A from Eqs.(33) and (34), we can find the rotational contributions 
to the well known Schwarzschild orbits for light and planets [36]. The geodesic 
Eqs.(33),(34) are valid only under the assumption of small test body approxi- 
mation (back reaction ignored, as stated in Sec. II, before Eq.(8)) 

Eq.(33) admits an immediate interpretation, albeit in optical parlance, as 
describing the motion of a particle in a "potential" (— l/2)n 2 v 2 and subjected 
to a "Coriolis" force 4j x (V x g), which would appear, for instance, in the 
absence of gravity (n = 1, h = 1) in a coordinate system rotating with angular 
velocity il = (1/2)(V x g). One might look at the first term in Eq.(33) as the 
"gravistatic" part and the second term as the "gravimagnetic" part. However, 
such gravimagnetic analogy practically works out only in the weak field limit of 
physically interesting metrics. The reason is that the coordinate transformation 
needed to achieve the spatially isotropic form in general entails use of implicit 
functions. To avoid this complication, one focuses attention only to large dis- 
tances on the equatorial plane (see Sec.V below). It would be interesting to 
compare and contrast the full optical analogy equations with the general rela- 
tivistic mechanical equations developed in Rcf.[39]. The gravimagnetic Eqs.(33), 
(34) do yield the known equations in the weak field, slow velocity limit [40] . In 
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CS modified gravity, the gravimagnetic analogy was first applied by Alexander 
and Yuncs [28] followed by Smith et al [29]. 

If one considers that gravity is given by the nonrelativistic, Newtonian law 
in an otherwise flat space (n — 1, h = 1) so that dA = dr, then one 



y2 _ 2M 



has from Eq.(33), the component of radial acceleration 

v x , 



dvj_ 

" - 9 + 2 



which leads to 



where 



A(f> SR = 



-9 



-.9 + 2 



v x f2 



pA 



1 d 



(v 2 ) 



M 



(35) 



(36) 



(37) 



2 dr r'- 

One can see how beautifully the effects of gravity and rotation, otherwise consid- 
ered as separate components in the literature [cases (ii) and (hi) in Sec. Ill], are 
synthesized into a single equation. We are now ready to apply Alsing Eqs.(33), 
(34) for calculating the phase shift in CS gravity. 

V. Phase shift in CS gravity 

To apply Eqs.(33), (34) we need to cast the KMT metric (6) in a spatially 
isotropic form. In general this is not possible, so for practical reasons, we 
restrict ourselves to equatorial plane (9 = tt/2) and to the weak field regime of 
the rotating source, r ^> 2M. Neglecting iv 2 , we have 



ds 2 



-h(r) 



dt + -rr^dcj) 



h(r) 
h(r) 

9<t> = ~- 



+ 1 



— dr 2 +r 2 d(f> 2 



1 



2M 
r 

u(r)r 2 

WT' 



(38) 
(39) 
(40) 



To achieve the isotropic form, we adopt isotropic radial variable p via $ = £ 
and ^7 = (1 + M/r)(dr/r), and obtain p = re~^~ , which yields r = p + M, 
to lowest order in — . Thus, in all, we have the expressions with n given by 
Eq.(32): 

M 2M 

$ p )~1-— ,h(p ~1 , 

p p 



2M 



2M\ 



n(p) ~ 1 + , ,g ~ - 1 + )(p+ M) 2 lu( p + M) 



P ) 

and the metric (38) for ^ <C 1 takes the form 



(41) 
(42) 



ds 2 = - 1 - 



2M 



[dt + g<pd(f>Y 



2M\ 
P 



)[dp 2 



(43) 
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With these at hand, and using Eq.(33), we can rewrite the radial and cross 
radial components of "acceleration" as follows: 



d 2 p 
dA 2 



dA 



dp 



4? 

dA 



(44) 



P dA 2 dAdA 



1 dp dg^ 
pdA dp ' 



(45) 



where (Ie^Vx ~~g) is the angular velocity of rotation of the gravitating mass 
M. The last term in Eq.(44) represents "Coriolis force". The first integral of 
Eq.(45) is 



d_ 

dA 



P 



dA 







dA 



9<p = 



(46) 



where t is an arbitrary constant of integration, which may be interpreted as the 
conservation of total "angular momentum" . Thus 



d0 
dA 



2-94 



(47) 



So far we have been talking in the language of optical analogue ("..") of 
mechanical quantities. Let us now translate them into mechanical quantities 
proper, quantities that arc actually measured in a gravity field. Thus we cal- 
culate the radial acceleration as follows. Converting A — > r by means of 
Eq.(34), we can rewrite Eq.(44) as 



dv± _ d 2 p ( d(f> 
dr dr 2 ^ \ dr 



(48) 



1 



d 

dp 



(49) 



where we have used in Eq.(44) the following: 



gether with 



dp 

717 



dr 
dA 



^ dr dA ) "™ 



to- 



2 V~h, being the angular velocity of rotation of M 



normal to the plane of the interferometer. We now need to calculate (3. Using 
Eqs.(29), (34) and Eq.(47) we find 



- d<t> d(j}dA 

(3 = 1 - 9 *Tt =1 - 9 *dAtt =1 - 9 * 



n 2p2 



(50) 



Since we disregard terms of the order of g^, being of the order of uj 2 , we end up 
with 

' r i + S ■ (51) 



r 1 



2 2 

n z p z 
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Then noting that v 2 = 2M/r ~ 2M/p, and using n from Eqs.(42) and (3 from 
Eq.(51), we can finally express v of Eq.(32) as 



v = 




n 2 p 2 



(52) 



On the other hand, we have the conserved proper energy and momentum for a 
test particle of rest mass m as follows [cf. Eq.(27)] 



E 



Ey/h 



(53) 



P = 



VT 



(54) 



such that the mass-shell condition may be written in a special relativistic form 

m 2 = hE 2 -<p- 2 p 2 = E 2 -p 2 . (55) 

(Note that, we could obtain, alternatively, the refractive index n by putting 
m = so that c = £ = c ^Vh = c n _1 .) Using Eq.(54), we can express v 2 in 
terms of p and $, and further introducing the quantum relation p = Hk, we get 
the curved space corrected Newtonian gravity term, A(j) g — > A0 g : 

' m 2 cl +p 2 $ 2 



~ 9 mp = - Ag 



= -9 



Am 2 

w 

Am 2 



1 



(56) 
(57) 



An interesting observation is that the relativistic second term is completely 
unaffected by gravity. Symbolically, we write the phase shift in CS gravity as 



CSG 



( —L 

{ dr 



1 

n 4 h 



pA 
Jw 2 



d_ 

dp 

= A^ + A^, 



2 2 

n v 



(n 2 Vhj {2Q n v) 



pA 
7m 2 



(58) 

(59) 
(60) 



where and A^> 2 are gravitational and Coriolis contributions respectively. 
This is the exact expression that contains terms of all orders in (i) and higher. 

All we have to do now is to collect the expressions for $, h, n, g v , v and ^M=- 
from Eqs. (41), (42), (49) and (52), insert them in Eq.(60), expand and find 
contributions to lowest order in (A). 



VI. Lowest order terms 
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To extract ostensibly first order terms, we ignore second and higher order 
terms in M, C\, C2 together with their various products. Then the shift A(f> 1 
to lowest order is as follows: 



A^ 



1 



h J \ Hv 2 J dp \ 2 



pA \ d ( n 2 v 2 



1 + 



,f4C 2 6C1 



I P 



+ 0(e 2 ) 



(61) 
(62) 



We can see how the term A^, which represents only A<p g in flat space, is 
modified by the quantities C\ , C2 and — . 

To calculate the Coriolis effect, we introduce the quantum relation E = Hgj 
so that 



hv 2 j 



(2Q n v) (^) = 2m n A<f>~ 1 . 



Thus we get from Eqs.(60) and (64) 



A0 2 



(n 2 y/hj 1 $- 1 (2Sa i ^) = + -y") (2u5f2 n A). 



Using the dispersion relation from Eq.(55) 



~2 w 



2-^2 



/i(p) ft.(p) I to 2 



(63) 
(64) 

(65) 
(66) 



in 2ujQ n A, we explicitly obtain the gravity corrected Coriolis term, A<fi c 



A0 C 



= 2ujfl n A 
M 

> I 1 + — 

P 



2AO r 



1 + t) 



Thus we have, from Eq.(65), 

A^ 2 



1 + 



AM 



Recombining it with Eq.(62), we get 



A0, 



CSG 



(A</> ff ) + 1 



4M 



(A0 C 



(67) 
(68) 



(69) 
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Recall that C\ has the dimension of length, C2 is dimcnsionless and I has the 
dimension of length, hence the CS effect £^j^- + i s dimcnsionless, as 

should be the case. Eq.(69) is the main result of our paper. 

As we see, the effect is coupled with I which is a sum of "angular momentum" 
P^^dk °f P ar ticles and of the source. The interferometer plane is usually placed 
at two vertical orientations so that the angular momentum of particles moving 
in the plane is practically zero. We are then left with only the source influence 
g^, which, when multiplied by C, becomes an effect in 0(e 2 ). Nevertheless, 
a number of observations can be made. It is immediate that in the limit of 
flat space (M = 0, C\, C2 = 0), we recover the individual special relativistic 
components Acf> g and A(f> cor derived by Anandan [14]. This means, in this 
limit, A<j) CSG — Ac/) SR . If C\, C2 = 0, but M 7^ 0, we recover the influence of 

Schwarzschild gravity embodied in the factors of , so that A(f> CSG — A(j) Sch . 
It is a new result by itself. If the black hole mass M — but C\, C2 7^ 0, we are 
left with the effect of a mere rotating coordinate system, which is of no physical 
interest. The remarkable result is that the parameters C\, C2 do not affect the 
Coriolis term at all, while it modifies only the pure gravity term A<f) g . This is 
because v, the carrier of C'i, C2 via g^, cancels out as shown in Eq.(64). 

Since the CS component of the phase shift depends crucially on the factor 
£, let us analyze what different values of it means. Although it is an optical 
constant, we see from Eq.(46) that it is also constant in proper time, 4^ = 0, 
because |j / 0. Note further that £ consists of two parts: A part p 2 |^, 
proportional to the usual particle orbital angular momentum, and a part g^ 
proportional to the angular momentum of the black hole itself. If i = 0, the CS 
effect vanishes altogether. This is the case when a radially freely falling particle 
acquires an angular velocity oc % near the source in the same sense as the 
rotation of the black hole, exactly similar to what happens in the Lense-Thirring 
drag [41]. On the other hand, if I = g^, then p 2 ^ = => & = 0, that is, the 
particle is at rest at constant (r, 9) relative to distant stars, and to remain at 
rest it has to swim against the spacetime drag. In this case, we see from Eq.(51) 
that there can be no 0(e) effect because the lowest order CS component would 
go like j j jXC~ 0(e 2 ), which we have agreed to ignore. Locally nonrotating 
[42] particles will have an angular velocity 4j£ = —gt4>/g<t><t> = but again 

the order of the effect would go like ui x C <~ 0(e 2 ). Thus the CS effect can 
appear only in 0(e 2 ). This is a bit unfortunate. The spin precession effect, as 
discussed by KMAT, seems better suited to study the CS effect. 

VII. Summary 

As shown by KMAT [33] , the constant parameter C2 in their solution is the 
observed (roughly constant) tangential velocity of particles moving in equatorial 
circular orbits in the galactic halo. Our aim here was to theoretically investigate 
the influence of such CS parameters on the quantum phase shift in thermal neu- 
tron interferometry. To achieve this goal, several steps were needed. We first 
established the stability of circular orbits in the CS theory, which is essential 
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for the physical viability of the model itself. The second step was to preview 
Anandan's arguments. In the third step, we familiarized the readers with the 
notion of curved spacetime optical-mechanical analogy, especially the reformu- 
lation of geodesic equations developed by Alsing [36] . (We relax the assumption 
that KMT solution (6) replaces that of Kerr but there are indications that the 
CS modification is a non-trivial perturbation of Kerr geometry [35]). In the 
final step, we applied these equations in the linear aproximation to obtain the 
effects on the phase shift contributed by inertial forces, special relativity, general 
relativity and the CS modified theory. 

Eq.(69) shows that, to leading order, the CS constants influence only the 
pure gravity term A</> g . However, the constant C\ is observationally clueless 
yet, but C2 already has a solid observable meaning. We can set C\ = without 
any loss of rigor, since u>(r) of Eq.(7) would still satisfy the core CS equation 
[r(r - 2M)w" + 4(r - 2M)J + 2uj = 0]. Stability of circular orbits is not 
hampered when C\ — 0. Hence, leaving aside relativistic terms 0(cq 2 ), we 

have A0 CSG ~ [1 + (a + /3) + af3] (A0 3 ) + (1 + 7 )(A0 cor .), where a = 

= ~ ' T = ' re P resen ting respectively the CS and Schwarzschild influences 
on the observed flat space effects of Newtonian gravity and inertial forces. 

One might argue that the mass M in question is actually that of the central 
galactic bulge and any possibility of interferometry is out of question. We wish 
to argue differently. Observationally, the frequency shifts in the HI radiation 
show that, in the halo region, C2 = v v /c is nearly constant at a value 7 x 10~ 4 
[43] . We may take this value of C2 as permanently fixed (up to observational un- 
certainty), like cosmological constant, and treat the KMT solution as completely 
determined like Schwarzschild-de Sitter solution. Then we can apply the solu- 
tion to study the nonlocal effect of this constant on the measurement of phase 
shifts for any rotating source, including Earth. To lowest order, and assuming 
w ~ -gfa in Eq.(42), we see that a = ^*£l = 4 Cf (1 + 2M) ~ 2 x 1CT 6 . The 

Schwarzschild (/3, 7) corrections arc of the order of — <~ 10~ 9 on Earth. Thus 

the current accuracy level measuring A</> g must be augmented by 10 6 times for 
CS effect and 10 9 times for the Schwarzschild effect. Whatever be the experi- 
mental scenario, our final result, Eq.(69), provides us with a clear view of the 
various effects within a single theoretical framework. It would be interesting to 
do analogous calculations for Dirac particles. But it would be highly nontrivial 
since the Dirac Lagrangian density would have to be coupled to the CS field, 
through the covariant derivative. Perhaps a starting point could be the analysis 
of Alexander and Yunes [31] but extending this to spinning particles would be 
very hard. 

Finally we would like to summarize the problematic characteristics of the 
KMT solution it has inherited from the non-dynamical framework. First, the 
solution is not a well accepted replacement of slow rotation limit of the Kerr 
solution. The KMT solution is a very special line clement that holds only for 
a particular choice of the CS coupling field. Unlike in the Kerr case, there is 
no Robinson's theorem [44] for rotating sources in CS gravity that could sug- 
gest that this solution is unique. Indeed, Yunes and Pretorius [35] have found 
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a family of new stationary axisymmetric solutions in the non-dynamical the- 
ory, which shows that solutions in this theory can not be unique. From this 
standpoint, the far-field solutions of Alexander and Yunes, and Smith et al are 
equally valid. Second, the KMT solution is incompatible with the dynamical 
formulation of CS modified gravity. This is because the energy contained in the 
KMT CS scalar is infinite. The Yunes-Pretorius solution bypasses the problem 
with the stress tensor and is compatible with dynamical theory. Third, the 
non-dynamical framework (upon which both the KMT solution [34] and the 
far-field solutions of Alexander and Yunes [27, 28] are based) is problematic for 
Schwarzschild black hole perturbation theory, as shown by Yunes and Sopuerta 
[32]. In this work, the authors showed that there cannot be single-parity per- 
turbations in non-dynamical CS modified gravity, which suggests that when a 
small object falls head-on into a supermassive black hole, the latter cannot be 
perturbed to linear order. This obstacle is partially resolved in the dynamical 
theory. Fourth, the non-dynamical framework is also problematic from the the- 
oretical standpoint. The main objection comes from the external, purely ad hoc 
prescription of the CS coupling field, which is not motivated by considerations 
from within the theory. Moreover, the non-dynamical framework also contains 
an additional constraint (the Pontryagin constraint), which suggests that the 
modified field equations are over-constrained. Indeed, this constraint affects 
the black hole perturbation theory (Yunes and Sopuerta [32]), and the search 
for exact solutions that could represent spinning black holes (Grumiller and 
Yunes [30] ) . If this is the case, then one might conclude that the non-dynamical 
framework can only have an ill-posed initial value problem [32]. 

Figure captions: 

Fig,l The height represents V" for small values of C\ ^ 0, C2 7^ along 
x and y axes respectively. The radius is arbitrary but fixed at any r = R. 
Negative values of V" indicate stability. 

Fig. 2. The de Broglie wavelets originating at A and interfering at B, one 
travelling along ADB and the other along ACB, the plane ADBC being rotated 
about the horizontal axis AO (Ref. [14]). 
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